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We theoretically study the DC Josephson effect of a semiconductor nanowire (NW) with strong spin-orbit 
interaction when a magnetic field is applied parallel to the NW. We adopt a model of single scatterer in a quasi- 
one-dimensional system for the case of short junctions where the size of normal region is much smaller than 
the coherent length. In the case of single conduction channel in the model, we obtain analytical expressions for 
the energy levels of Andreev bound states, E„, and supercuiTent /, as a function of phase difference ip between 
two superconductors. We show the Q-n transition by tuning the magnetic field. In the case of more than one 
conduction channel, we find that E„{—ip) ?t EiXf) by ths interplay between the spin-orbit interaction and Zeeman 
effect, which results in finite supercurrent at i/j = (anomalous Josephson current) and direction-dependent 
critical current. 

PACS numbers: 
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For spin-based electronics, spintronics, the manipulation of 
electron spins in semiconductors is an important issue.' The 
strong spin-orbit (SO) interaction in narrow-gap semiconduc- 
tors, such as InAs and InSb, has attracted a lot of interest 
in this context.^ Nanowires (NWs) of such materials have a 
great potential for the application to the spintronic devices 
and also to quantum information processing by utilizing the 
SO interaction.^^ Indeed, the electrical manipulation of sin- 
gle electron spins was reported for quantum dots fabricated 
on the NWs.''"^ Recently, the proximity effect was intensively 
examined when NWs are put on superconductors, for the 
search of Majorana fermions by the combination of SO inter- 
action and magnetic field. ^ The DC Josephson effect was also 
studied when the NWs are connected to two superconductors 
(S/NW/S junctions). "^'^ In this paper, we theoretically study 
the effect of strong SO interaction on the Josephson current 
through the NWs in a magnetic field, which is closely related 
to the recent experimental results. — 

The Josephson current through mesoscopic systems of nor- 
mal metal or semiconductor has been studied for a long time. 
At the interfaces between the normal systems and supercon- 
ductors, the Andreev reflection takes place in which an elec- 
tron (a hole) is conveited to a hole (an electron). As a re- 
sult, the electron and hole are coherently coupled to each 
other, forming the Andreev bound states in the normal re- 
gion around the Fermi level within the superconducting en- 
ergy gap Aq. They have discrete energy levels £„ (Andreev 
Ievels)ji^4i^ In the presence of phase difference (f between the 
superconductors, the supercurrent Kip) is carried by the An- 
dreev bound states when the length of the normal region L is 
much smaller than the coherent length ^ (short junction) ji^— 
^ - hvp/iTrAo) = ^0 for ballistic systems and ^ - (^oO'^^ 
for diffusive ones, where vp is the Fermi velocity and / is the 
mean free path. The supercurrent is simply written in terms 
of the transmission probability r„ for conduction channel n 
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In Josephson junctions of superconductor / ferromagnet / 
superconductor (S/F/S), the 0-n transition was observed and 
intensively studiedi2£i2i In the 7r-state, the free energy is min- 
imal at (fi = TT, which stems from the Zeeman splitting by 
exchange interaction in the ferromagnet. The splitting makes 
the spin-dependent phase shift for electrons in the propaga- 
tion through the ferromagnet. Since the Andreev bound states 
consist of an electron with spin cr and a hole with spin -cr, the 
Andreev levels are dependent on the spin in the feiTomagnet. 
The Q-7T transition was observed when its thickness is grad- 
ually changed, as a cusp of critical current^^ A similar tran- 
sition was recently observed in S/NW/S junctions with fixed 
length when the Zeeman splitting is tuned by applying a mag- 
netic field parallel to flie NWM 

The effect of SO interaction is another interesting subject 
for the DC Josephson effect. It was investigated by a lot of 
theoretical groups, for normal metal with magnetic impuri- 
ties,^^ two-dimensional electron gas (2DEG) in semiconduc- 
tor heterostructures,^^^^ open quantum dots (QDs),^"* QDs 
with tunnel baiiiers,''''^*' carbon nanotubes,^' quantum wires 
or NWs,"*^^^ and others.'*'' Even in the absence of magnetic 
field, the SO interaction splits the spin-degeneracy of the An- 
dreev levels when the phase difference if is finite.^** "^^ In the 
short junctions, however, the splitting is not observed unless a 
weak energy-dependence of the scattering by the SO interac- 
tion is taken into account. In this case, the supercurrent is 
given by eq. ([Til, irrespectively of the SO inteaction. 

The coexistence of the SO interaction and Zeeman eff'ect 
induces a supercurrent at ^ = 0, so-called anomalous Joseph- 



son current 
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The anomalous cuiTent flows in the 



i/jQ-state in which the free energy has a minimum at (/? = (^o 
{i^ 0, tt)."*^ The anomalous Josephson current was predicted 
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when the length of normal region L is longer than or compa- 
rable to the coherent length ^. Krive et al. derived the anoma- 
lous current for long junctions (L » ^) with a single conduc- 
tion channel. Reynoso et al. found the anomalous current 
through a quantum point contact in the 2DEG for L > 
They also showed the direction-dependence of critical current 
when a few conduction channels take part in the transport. 
In the experiment on the S/NW/S junctions,'^ the direction- 
dependent supercurrent was observed for samples of L > ^ 
in a parallel magnetic field, besides the above-mentioned O-tt 
transition. This should be ascribable to the strong SO inter- 
action in the NWs although the anomalous Josephson current 
was not examined. 

In this paper, we study the properties of the supercurrent 
through NWs with strong SO interaction, focusing on the case 
of short junctions. We elucidate the anomalous Josephson cur- 
rent and direction-dependent critical current, based on a sim- 
ple model. In our model, both elastic scatterings by impurities 
and strong SO interaction in the NWs are represented by a sin- 
gle scatterer in a quasi-one-dimensional system. The number 
of conduction channels is unity or two. The Zeeman effect 
is taken into account by the spin-dependent phase shift in the 
propagation through the system. 

First, we analyze the model with N - I. We calculate the 
Andreev levels E„ as a function of ip, which yields the super- 
current I{ip) via eq. (|5]l. We obtain an analytical expression 
for I{ip) in the absence of SO interaction and clearly show the 
O-TT transition when the magnetic field is tuned. The SO inter- 
action does not change the supercurrent qualitatively in this 
case. We still find the relation of E„(-(p) - E„(ip) and that the 
free energy is minimal at = or tt. We observe no anoma- 
lous supercuiTent. 
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FIG. 1 : Our model for a semiconductor nanowire (NW) connected 
to two superconductors. The NW is represented by a quasi-one- 
dimensional system along the x direction, with a single scatterer at 
X = Xq to describe the impurity scattering and strong SO interaction. 
The Zeeman effect is taken into account by the spin-dependent phase 
shift for a magnetic field applied in the x direction, (a) Schematic 
view of the model. The pair potential in the NW is Aoe"''- at x < 
and Aoe'^"* at L < x by the proximity effect, and at < x < L. 
(b) The scatterer at x = xo is represented by the scattering matrix 5 . 
ae(h) and 6e(h) denote incoming and outgoing electrons (holes) with 
respect to the scatterer, respectively, whereas b'^^^^^ and a'^^^^ are those 
with respect to the boundary at x = and L. Ep is the Fermi energy. 



Next, we examine the model with = 2, in which the inter- 
channel scattering takes place at the scatterer The scattering 
is represented by a random matrix of the orthogonal ensemble 
in the absence of SO interaction and that of the symplectic en- 
semble in the strong limit of SO interaction. The ensembles 
are interpolated for the intermediate strength of SO interac- 
tion. We show that the interplay between the SO interaction 
and magnetic field results in the breaking of E„{-(p) - E„{ip) 
for the Andreev levels, in contrast to the case of N - 1 . This 
leads to the anomalous Josephson current and the direction- 
dependent current. We believe that our model should be use- 
ful to understand the experimental results on S/NW/S Joseph- 
son junctions in which a few conduction channels exist in the 
NWs,'"* as an approach from the limit of short junctions. 

The organization of this paper is as follows. In Sec. 2, we 
explain our model of single scatterer for the S/NW/S Joseph- 
son junctions and calculation method of the Andreev levels 
and supercurrent. Analytical results are given in Sec. 3 for 
the case of single channel, whereas numerical results are pre- 
sented in Sec. 4 for the case of two channels. The last section 
(Sec. 5) is devoted to the conclusions and discussion. 



II. MODEL AND CALCULATION METHOD 

In this section, we explain a model of single scatterer in a 
quasi-one-dimensional system for the NWs. The Bogoliubov- 
de Gennes (BdG) equation is introduced for the calculation 
of the Andreev bound states. The formulation of solving the 
BdG equation is given in terms of the scattering matrix. 



A. Model and BdG equation 

Figure [Tfa) shows our model in which a NW along the x 
direction is connected to two superconductors. At x < and 
L < X, the Cooper pairs penetrate into the NW by the prox- 
imity effect. They are transported through the normal region 
of < jc < L. The Zeeman effect is taken into account for a 
parallel magnetic field B which is so weak as not to break 
the superconductivity. The Hamiltonian is H = Hq + Hz, 
with Hq - 12m + Vconi(y,z) + Vscatt and Zeeman effect 
Hz - gfiB^ ■ (t/2, using effective mass m, ^-factor ^ (< 
for InAs and InSb), Bohr magneton /yb, and Pauli matrices &. 
We take the spin quantization axis along the magnetic field (x 
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direction). Vconf describes the confining potential in y and z 
directions. The number of conduction channels is = 1 or 
2. Vscatt represents both the elastic scattering and SO interac- 
tion at a single scatterer at x - xq. Its explicit form is given 
by the scattering matrix in the next subsection. We assume 
that L <K ^0 and there is no potential barrier at the bound- 
aries between the normal and superconducting regions. Both 
the Zeeman energy Ez - \gfJ.BB\/2 and the superconducting 
energy gap Aq are much smaller than the Fermi energy Ep. 

To obtain the Andreev bound states, we solve the BdG 
equation/''''*** 
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where - (i/'e+, lAe-)^ and i^h = ('Ah+>'Ah-)^ are the spinors 
for electron and hole, respectively. The energy E is measured 
from the Fermi level Ep. A is the pair potential in the spinor 
space 
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where g - -/o",,.'*^ We assume that A(x) is Aoe"^'- at < 0, 
at < ;c < L, and Aoe"^" at L < x. The BdG equation 
determines the Andreev levels E„ {\E„\ < Aq), as a function of 
phase difference between the superconductors ip - ipi^- (/jr. 

If the BdG equation has an eigenenergy E„ with eigenvec- 
tor (l^e.n, <Ah,«)^i it also has the eigenenergy -E„ with eigen- 
vector (^^'hn'^^en)^- short junctions of L <k ^o, the number 
of Andreev levels is given by AN; 2N positive levels and 2N 
negative ones when the number of channels is (2N if the 
spin degrees of freedom is included). The ground state energy 
is given by 
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where the summation is taken over all the positive Andreev 
levels, E„{(p) > 0. The contribution from continuous levels 
(\E\ > Ao) is disregarded in eq. (|4]l, which are independent of tp 
in the short junctions.'* At zero temperature, the supercurrent 
is calculated as 



2e dEg^ 
fi dip 



e ' dEn 
dip 
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The current is a periodic function for -n < ip < n. The max- 
imum (absolute value of minimum) of liip) yields the critical 
current /c_+ (/c,-) in the positive (negative) direction. 

The symmetry of the BdG equation should be noted here. 
We denote the matrix on the left side of eq. ^ by 'H(ip). 
In the absence of Zeeman effect, T'J-(((p)'J ' = 'J-({-(p) 
with the time-reversal operator T = -i&yK for spin- 1/2 
particles. If "/-{((p) has an eigenenergy E„ with eigenvector 
ii^e,„,i^h.n)^ , 'J^i-f) has an eigenenergy E„ with eigenvec- 
tor Tiip^ n, if^h,nV- Thus the Andreev levels satisfy the rela- 
tion of E„((p) = E„{-(p). In the absence of SO interaction, 
K'HiipW-^ = 'H(-ip). Then we derive that E„((p) = E„i-ip) 
in the same way. The relation does not always hold in the 
presence of both SO interaction and magnetic field. 



B. Formulation using scattering matrix 

The BdG equation in eq. ^ is written in terms of the scat- 
tering matrix.'* 

First, we represent the effect of Vscatt by 5p (p = e, h): At 
the scatterer at x = xq, an electron is scattered to an electron 
by 5e and a hole is scattered to a hole by Sh- They connect 
the amplitudes of incoming particles of conduction chan- 
nels with spin ±, (flpL, dpRp', and those of outgoing particles. 



«pL 
«pR 



(6) 



Se and Sh are 4N x 4-N matrices and related to each other by 
Se(E) = S^{-E). On the assumption that they are independent 
of energy E for \E\ < Aq and thus 5e = 5^, we denote = S 
and Sh = S*. S is conventionally written by reflection and 
transmission matrices: 



S = 



(7) 



i, S satisfies that 



In addition to the unitarity, S' S 

g^'^Lg, = g^f^Rg, and t^^ - g^ ii^Rg when the time reversal 
symmetry holds. 

Second, we describe the transport of an electron (a hole) in 
< X < xo and xq < x < L by scattering matrix tb (t^), 
considering the Zeeman effect Hz- For \E\ Ep, we use a 
linearized dispersion relation, E = +hvp(k - kp) for k > Q and 
E - -fivp{k+kp) for ^ < 0, where kp is the Fermi wavenumber 
For spin cr = +1, the wavefunction is i^f^ h ^ where k = 
kp + (E + Ez)l(fivp) fork>Q and k ^ -kp - {E + Ez)l{tivp) 
for k < O."*^ In the Andreev bound states, a pair of right-going 
(left-going) electron with spin cr - +\ and left-going (right- 
going) hole with spin cr -+\ acquire the phase of ifl^L with 



\gWB 
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Y"''- 2nvp 



■xo 



in the propagation at < x < xq, and ±0br with 



1 



(L - xo) 



(8) 



(9) 



in the propagation at xq < x < L. We can safely disregard 
the phases of 2Exo/(hvp) and 2E(L - x^))/{fivp) since l^l < Aq 
and L <k ^o- In the case of two channels (A^ = 2), the Fermi 
velocity is different for the channels. We neglect the differ- 
ence because it would not qualitatively change our numerical 
results in Sec. 4. 

The spin-dependent phase shift is represented by the scat- 
tering matrix for an electron 



= Tb 



where (b'^i^, b'^^)^ are amplitudes of outgoing electrons 
Oorx = L[Fig.[nb)]. It is 



Tb = 
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where 



III. CALCULATED RESULTS FOR SINGLE CHANNEL 
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a/2 
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and i is the N x N unit matrix. In our model, the Zeeman 
effect is characterized by two parameters. One is its strength. 
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where E-rh - fivp/L is the Thouless energy for the ballistic 
systems. The other is an asymmetry between Obl and Obr, 
Ob - Obl/Obr - xo/(L - xq). We fix at = V2 for the 
calculations in this paper 

Third, the Andreev reflection at x = and L is described by 
scattering matrix rhe for the conversion from electron to hole 
and Teh for that from hole to electron. When an electron with 
spin cr is reflected into a hole with -cr, it is written as 
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with a - arccos(£'/Ao). When a hole is reflected to an elec- 
tron, it is 



'•eh I 5^ 



with 



'"eh = 1 
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We assume that the channel is conserved at the Andreev re- 
flection in the case of N = 2. The normal reflection can be 
neglected in our case without potential barriers at the bound- 
aries of X = and 

The total scattering matrix is obtained by the product of S , 
tb, he, and reh- The BdG equation yields'^ 



det(i - fBrehfJj5*f^rhetij5) = 0, 



(18) 



which determines the Andreev levels E„((p). In the absence of 
magnetic field, eq. (fTST l is simply reduced to - 



det 



1 -1^1 -fiprLRsin 
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= 0. 
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In this case, the Andreev levels are represented by the trans- 
mission eigenvalues of f^^ii^R- They are two-fold degenerate 
reflecting the Kramers' degeneracy. Thus the Andreev lev- 
els E„((p) are not split by finite (p in spite of the broken time 
reversal symmetry. 



In this section, we present the analytical results for the case 
of single channel (A^ = 1). In this case, the reflection matrices 
in eq. (|7]i are generally written as tl = VT-T and tr = 
- Vl - r where T is the transmission probability of the 
scatterer The phase factor e'^ is cancelled out in eq. ( fTSl l. The 
transmission matrices are 



fLR 



p'lSO 



-"/so 
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and 



where 



fRL 



u. 



cos(0/2) sin(6l/2) 
-sin(6'/2) cos(6i/2) 



(20) 



(21) 



(22) 



The SO interaction is described by the rotation around an axis 
(effective magnetic field) tilted from the spin quantization axis 
(x direction) by 6. 7750 characterizes the strength of SO inter- 
action; < //so ^ The matrix S belongs to the orthogonal 
ensemble of random matrix theory for //so = (no SO inter- 
action) and to the symplectic ensemble for t/so = (strong 
limit of SO interaction). 



A. In absence of SO interaction 

We begin with the case without SO interaction (t/so = 0). 
The Q-n transition is elucidated as a function of magnetic field. 

By solving eq. ( fTSl ), we obtain the analytical expression for 
the Andreev levels 



E^±(tp) 



Ei±{(p) 
Ao 



J H- arccos |^± + 6b + T cos </3)/2j j , 



(23) 



> - arccos |^+ + 6b + T cos vp)/2j j , 



(24) 

where 6b - {I - T)cos[6B{aB - 1)/(Q'b + 1)]- The subscript t 
(i) indicates the state of electron spin 0- - +1 {(t = -I) and 
hole spin cr - -I (cr = +1). The sign + in eqs. (|23T l and ( l24l i 
corresponds to the positive or negative energy at 6b - 0, i.e. 
Ei+((p) - Ei+{ip) - -Ei^{(f) - -Ei^iif) because of the spin 
degeneracy in the absence of magnetic field. As mentioned in 
Sec. 2.1, positive and negative levels appear in pairs; E-^+iip) - 
-Ei-((p) and Ei+{if) - -E-i-{ip) even for 6b + 0. The Andreev 
levels are an even function of ip, En{-ip) - E„{(f) for n = (T 
±), (i ±). 

We plot the Andreev levels in eqs. ( |23T l and ( l24l i as a func- 
tion of If in Fig. |3a). The magnetic field gradually increases 
from = to 71. We find three regimes for 6b- 

(I) In the regime of < 6b < 0^*, E^+,Ei+ > and 
Ei-,Ei^ < 0. A weak magnetic field splits the Andreev levels 
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for spin t and J,. With an increase in 9b, the sphtting increases 
and finally E^^ - Ey^ - si ip - ±n when Og - 0^*. The 
ground state energy E^s in eq. (|4| takes a minimum at = 
(0-state) in this regime. 



(II) When < 6b < 0^^', E^+itp) and Ei^{tp) intersect at 
(fi - +(p[ and £■ = 0; Ei+{(pi) - Ei^{(pi) - 0. tpi satisfies the 
condition of 



0(2) 



T cos (fii + cos 0B + Sb(Ob) - 0. 



(25) 



With increasing Ob, the intersections move from ±n (6b = 0g ') 
to {9b = 0^'). 0^'' and (f^^ are determined from cos 9b + Sb - 
T and cos 9b + 6b - -T , respectively. 

(Ill) When 9f < 9b < n, £'i+,£'i- > and £'t+,£'t- < 0. 
In this regime, Egs is minimal at = tt (;r-state). The Q-n 
transition takes place at a critical value of 9b in regime (II). 

The behavior of E„{(p) from 9b - ti Xo In is similar to that 
from 9b = tt to in Fig.|2ja). (They are precisely identical to 
each other in the case of xq = L/2.) 

The supercurrent I{ip) is evaluated using eq. (2). In regime 
(I), it is given by 




n<f) = — cos 
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T sin If 
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This equation coincides with eq. ([T]l in the absence of mag- 
netic field {9b = 0). I{(p) oc sin (fi for T «; 1, which is typ- 
ical for the 0-state. In regime (II), I((p) is discontinuous at 
(fi - +(fi . It is written as 
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FIG. 2: Calculated result for the model of single conduction channel 
in the absence of SO interaction (A' = 1, eso = 0). The transmis- 
sion probability at the scatterer is T = 0.8. (a) Andreev levels as a 
{\(p\ < (fii ) function of the phase difference (p between the superconductors, Ef+ 

(27) (solid lines) and (broken lines). The magnetic field is 6^ = 
(19^1 > Vi) (left upper), 0. l/r (left middle), OAk (left bottom), OJn (right upper), 

0.9n (right middle), and n (right bottom). At 6b = 0, solid and broken 



In regime (III), the supercurrent is continuous and given by 
eAo . 19b\ T sinip 



eAo . (9b\ 

sm — 

n \2I 



^J2 - 2{6b + T cos ip) 



(28) 



I((p) oc - sin(/3 for T <K 1, which is a character of the 7r-state. 
Figure |2jb) exhibits the supercurrent I{(p) as a function of (p. 
Note that the supercurrent satisfies /(-^) = -lif) because 
Egsi-if) = Eg,((p). 

In Fig. |2fc), we plot the critical current as a function of 
magnetic field. The critical current /c,+ in the positive di- 
rection is identical to 1^- in the negative direction. With in- 
creasing 9b, /c,± decreases first and turns to increase showing 
a cusp at the critical point of O-tt transition. The critical point 
is around 9b ~ 7t/2, or Ez ~ Ej^ from eq. ( fT3] l. 



B. In presence of SO interaction 

In the presence of SO interaction (?7so ^ 0), the transmis- 
sion matrix in eq. (l20l l is rewritten as 



lines are overlapped to each other; Ef+ = At dg - n, = E-^_ 
and = Ei-. (b) Supercurrent through the NW. The magnetic 
field is 6b = (solid line), QAn (broken line), Q.ln (dotted line), and 
n (dash-dot-line), (c) Critical current 7^.+ as a function of magnetic 
field, Qg. The current in the positive direction 1^ + is identical to that 
in the negative direction 



where eso 



sin^(77so) sin^ 9 and = 
Thus the spin-flip proba- 



fLR 



(29) 



arccos [cos(/7so)/ VI - eso]- 
bility is equal to eso- If the effective magnetic field of SO 
interaction is parallel to the magnetic field {9 - 0), eso = 
and no spin flip takes place. Since the phase is irrelevant 
to the calculation in eq. ( fTSl ), the effect of SO interaction is 
described by single parameter eso. 

We calculate the Andreev levels E„{ip) by solving eq. ( fTSl) . 
Figure[3la) shows the levels as a function of i/j when eso = 0.2. 
The magnetic field gradually increases from 0^ = to n. Since 
the spin states are mixed by the SO interaction, t and J, are 
not good quantum numbers. However, the influence of the SO 
interaction is inconspicuous in the case of single channel. 

In the absence of magnetic field {9b - 0), all £,,(1^) are in- 
dentical to those in Fig.|2a) with eso - because the eigen- 
values of ^LR^LR are T in eq. ( fT9] l, irrespectively of eso, when 
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In the case of single channel, the relation of E„(-ip) - E„{(p) 
holds even in the presence of SO interaction. In consequence 
the supercurrent satisfies I{-(f) = -lif) in Fig. IHb). We 
do not observe an anomalous Josephson current or direction- 
dependent critical current: I{if = 0) = and /c_+ - 1^-. 



IV. CALCULATED RESULTS FOR TWO CHANNELS 

In this section, we examine the case of two conduction 
channels in a NW (A^ = 2) by numerical calculations. 
We show an anomalous Josephson current and direction- 
dependent critical current, in contrast to the case of single 
channel. 

For the scattering matrix S in eq. dTji for a single scatterer at 
X - xo, we prepare random samples following the orthogonal 
ensemble in the absence of SO interaction and the symplectic 
ensemble in the strong limit of SO interaction. For the inter- 
mediate strength of SO interaction, the ensembles are interpo- 
lated with a parameter /?so (0 < pso ^1), using the method 
given in Appendix, pso - for the orthogonal ensemble and 
Pso = 1 for the symplectic ensemble. 

First, we present the calculated results for a sample in the 
absence of SO interaction (pso - 0) and that in its presence 
ipso + 0). Then the random average is taken for the latter 
case. 



FIG. 3: Calculated result for the model of single conduction channel 
in the presence of SO interaction {N = \, 6so = 0.2). The transmis- 
sion probability at the scatterer is T = 0.8. (a) Andreev levels E„ 
as a function of the phase difference if between the superconductors. 
The magnetic field is 6*8 = (left upper), 0. l;r (left middle), 0.4;r 
(left bottom), O.T/r (right upper), 0.9w (right middle), and n (right 
bottom). At ds = 0, two lines are overlapped to each other, reflect- 
ing the Kramers' degeneracy, (b) Supercurrent /((/)) through the NW. 
The magnetic field is flg = (solid line), 0.4;r (broken line), 0.7;r 
(dotted line), and n (dash-dot-line), (c) Critical current 4 + as a func- 
tion of magnetic field, ffg. The current in the positive direction /c,+ is 
identical to that in the negative direction /(.,_. 

fLR is given by eq. ( |29] l. Thus the Andreev levels are not af- 
fected by eso ■ 

The magnetic field splits the Andreev levels. The behavior 
of the levels with Qb is similar to that in the case of eso - 0. In 
regime (I), the splitting is enlarged with an increase in Qb- The 
0-state is realized here. In regime (II), the crossing points ±ip\ 
of two levels at £ = move from +-n toward 0. The equation 
for ^1 is modified to 

Tcos^Ji +cos0b + (Jb(0b) = -TTeso sin | J sin | "^^^ | 

\aB + \] yo-B + l/ 
(30) 

by the SO interaction. 6*^^^ and Q^'^ satisfy eq. (O at ip\ = ±ti 
and ip\ - 0, respectively. We do not observe regime (III) in 
Fig.|3la) since there is no solution for 4^ in this case. 

The supercurrent is shown in Fig. |3lb) as a function of 
whereas the critical current is in Fig. |3jc) as a function of Qb- 
They are qualitatively the same as those in the case of eso - 0. 



A. In absence of SO interaction 

Figure |4ja) shows the Andreev levels as a function of ip, 
for a sample in the absence of SO interaction {pso - 0). The 
magnetic field gradually increases from Qb = to tt. The 
levels are labelled as E^±i (solid line) or Ei±i (broken line) 
with / = 1,2 in the same manner as in Sec. 3.1. They appear 
in pairs, E^^^ii^p) = -Ei^ii^p) and ^x+j^^p) = -E^^-iiip). 

When Qb - 0, the Andreev levels are spin-degenerate, 
E^±iW) - Ei±i((p). When Qb + 0, we find three regimes for Qb 
just as before. In regime (I), £'i-+/, E^^i > and /i-j--/, Ei^i < 
although the levels are split by the Zeeman eff'ect. The ground 
state energy has a minumum at(f = Q (0-state). In regime (II), 
intersections between Ei+i{(p) and £'|_,(i/3) appear at £ = for 
/ = 1 only, or both of / = 1,2. In regime (III), £'j±,((^) > and 
Ei±iitp) < 0. The TT-state is realized in this regime. 

In Fig. 1Mb), the supercurrent I(ip) behaves almost in the 
same way as in Fig. 2(b) though I{(f) shows four disconti- 
nuities in the case of four intersections in regime (II). Figure 
|4fc) shows the critical current /c,±, which displays cusps corre- 
sponding to the 0-;7r transition at a critical values of Qb- The re- 
lation of E„((p) = E„(-ip) holds, which yields I{ip) - -I{-(f). 
Thus we do not observe the anomalos Josephson current or 
direction-dependence of /c,±. 



B. In presence of SO interaction 

In Fig. 13 a), we present the Andreev levels for a sample in 
the presence of SO interaction ipso = 0.3). The interplay 



7 








5 » ^ ^ 









I 

■^-1 



/ ■■•\ • 




— 1 • J 

\ 1 w 


\ y 








n -n 




(b) 



1 
-1 



1 \ / 
\ * x 

v ^ 


/a ^ 1 

t. "J 

''<:■:/ 


n C 





(p 




FIG. 4: Calculated result for the model of two conduction channels 
in the absence of SO interaction (A' = 2, pso = 0). A sample is ex- 
amined for the scattering matrix 5 at the scatterer. (a) Andreev levels 
as a function of the phase difference between the superconductors, 
E-^+i (solid lines) and (broken lines) with i = 1,2. The magnetic 
field is (Js = (left upper), 0. l;r (left middle), QAn (left bottom), Q.6jt 
(right upper), 0.8;r (right middle), and n (right bottom). At Gg = 0, 
solid and broken lines are overlapped to each other; E'j±j = -Ej^±,. At 
6b = ?r, E^+i = Ef-i and E^+j = E^-j. (b) Supercurrent I((p) through 
the NW. The magnetic field is fg = (solid line), OAn (broken line), 
0.8;r (dotted line), and n (dash-dot-line), (c) Critical current 1^. + as a 
function of magnetic field, 9b. The current in the positive direction 
Ic,+ is identical to that in the negative direction . 



FIG. 5: Calculated result for the model of two conduction channels 
in the presence of SO interaction {N = 2, pso = 0-3). A sample is 
examined for the scattering matrix 5 at the scatterer. (a) Andreev 
levels E„ as a function of the phase difference between the super- 
conductors. The magnetic field is = (left upper), 0.1:rr (left mid- 
dle), OAn (left bottom), 0.6n (right upper), 0.8;r (right middle), and 
n (right bottom). At Ob = 0, two lines are overlapped to each other, 
reflecting the Kramers' degeneracy, (b) Supercurrent I(ip) through 
the NW. The magnetic field is Sg = (solid line), OAn (broken line), 
O.Stt (dotted line), and n (dash-dot-line), (c) Critical current as a 
function of magnetic field 9b, 4.+ in the positive direction (solid line) 
and in the negative direction (broken line). 



between the SO interaction and Zeeman effect leads to a qual- 
itatively new situation. 

In the absence of magnetic field {6b = 0), the Andreev lev- 
els E„((p) are two-fold degenerate. The Kramers' degeneracy 
is not removed by finite ip, as discussed at the end of Sec. 2.2. 
In the presence of magnetic field, they are split and show that 
E„((p) + E„{—(f). The SO interaction mixes different conduc- 
tion channels in a spin-dependent way to form the Andreev 
bound states. This breaks the relation of E„{ip) - E„{-(p) 
when 6b + 0. Roughly speaking, we can identify three 
regimes for < flg < ;7r, as in Fig. H) The 0-state appears 
in regime (I) at 6b ~ 0, whereas ;r-state is realized in regime 
(III) at Qb ~ TT. The Q-n transition seems to take place in the 
intermediate regime. However, we do not observe intersec- 
tions between the Andreev levels at £ = in regime (II). This 
is due to the anti-crossing of the levels by the SO interaction. 

In Fig. |3b), we show the supercurrent I{(p) as a function 



of ip. I((p) oc sin ip nt dB = and - sin ip at 6b - n, which 
are typical behaviors in the 0-state and 7r-state, respectively. 
We do not observe the discontinuity of l(ip) in regime (II), 
reflecting the absence of intersections of Andreev levels. 

It should be stressed that / at = in Fig. [Sfb), indi- 
cating the anomalous Josephson current. We plot the anoma- 
lous supercurrent 7(0) as a function of 6b in Fig. |6] Since 
Essif) + E„si-ip) for the ground state energy, I(p) + -l(-p) 
for the supercurrent. This results in finite 7(0). In other words, 
Egfi((p) takes a minumum at (po 0) in regime (I) except 
6b - 0. Thus the so-called i^o-state is realized, where 7(0) 
in eq. (|5]l. Similarly, E^sitp) is minimal at ^) in regime 

(III). Our numerical result indicates a discontinuous change 
of (/)() at a value of 6b in regime (II), from ipo ~ to (po ~ n. 
Therefore, there is a well-defined transition from 0-like-state 
to ;r-like-state at the critical value of 6b- 

In addition to the anomalous Josephson current, the critical 
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FIG. 6: Supercurrent at <p = 0, 1(0), as a function of magnetic field 
dg, for tfie model of two conduction channels in the presence of SO 
interaction (N = 2, pso = 0.3). The sample is the same as that used 
for Fig. |5] 
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current depends on its direction. Figure |3c) shows the current 
/c + in the postive direction (solid Hne) and 1^- in the negative 
direction (broken line), as a function of 9b- Both 1^,+ and /c__ 
show sharp changes from a decreasing function to an increas- 
ing one around 9b = 7!'/2 and 37t/2. The values of 9b at the 
cusps are different for /c,+ and 1^-, which are located below 
and above the critical value of the above-mentioned transition, 
respectively. 

In Figs. |5] and |6] we have presented the results for a sample 
for the scattering matrix S at the scatterer when pso = 0.3. 
We perform the numerical calculations for 100 samples and 
take a random average over the samples. 

Regarding the anomalous Josephson current, 7(0), we plot 
</(0)) and V<[A/(0)]2), where A/(0) = /(0)-</(0)), as a func- 
tion of 9b, in Fig.|2ja). pso = 0.3. The average of anoma- 
lous current, {/(O)), is almost zero since it is positive or neg- 
ative depending on the samples. Its fluctuation ^J{[AI(0)]^} 
yields an estimated anomalous current, which is of the or- 
der of O.leAo/^. It is zero in the absence of magnetic field 
{9b - 0), increases with 9b, and becomes maximal around 
9b = TT. Then it decreases with 9b until 9b ~ 27r. 

Figure Hb) shows ^J{[AI(0)]^} at 9b ^ tt (solid line) and 
Ob - O.Stt (broken line), as a function of the strength of SO 
interaction, pso- The anomalous Josephson current increases 
almost linearly with pso for small pso- It should be observ- 
able when Pso ^ 0.05. 

Next, we examine the direction-dependence of the critical 
cuiTent, 6Ic = Ic,+ - Ic,-- Figure |3c) plots its random average 
{6Ic) and fluctuation ^J{[A(6Ic)V}, where A(6Ic) = Sic - (Sic), 
pso = 0.3. The random average yields (i5/c) ~ 0, whereas its 
fluctuation is of the order of O.leAo/ft at n/2 < 0^ < 37r/2, 
where the 7r-like-state appears. 

In Fig. Etd), we show ^/{[A(5IJp} at 9b = n (solid line) 
and 9b = O.Stt (broken line), as a function of the strength 
of SO interaction, pso- The direction-dependent supercurrent 
could be observed when pso ^ 0.05. 



FIG. 7: Calculated result for the model of two conduction chan- 
nels (A' = 2). The random average is taken for 100 samples for 
each strength of SO interaction pso- (a) Average of the supercur- 
rent at (/) = 0, {I{(fi = 0)), as a function of magnetic field 9^. Er- 
ror bars represent the average of the fluctuation, -^([A/(0)]^), where 
A/(0) = /(O) - </(0)). Pso = 0.3. (b) V<[A/(0)]2) at = ;r (solid 
line) and 6b = O.Sn (broken line), as a function of the strength of 
SO interaction, psQ. (c) Average of the direction-dependence of the 
critical cuiTent, SI^ = /c,+ -/c-, as a function of magnetic field 9b. Er- 
ror bars represent the average of the fluctuation, -\/({K(6I^yp}, where 
A((5/,) = (5/, - (Sh). Pso = 0.3. (d) ^J{[A(Shm at 6»s = tt (solid 
line) and 9b = O.Sn (broken line), as a function of the strength of SO 
interaction, pso- 



V. CONCLUSIONS AND DISCUSSIONS 

We have studied the DC Josephson effect in S/NW/S junc- 
tions in the presence of strong SO interaction in the NWs and 
Zeeman effect in a parallel mangetic field. We have examined 
a simple model of single scatterer in a quasi-one-dimensional 
system for short junctions where the length of the normal re- 
gion is much smaller than the coherent length (L <K ^). For the 
case of single conduction channel, we have obtained analyti- 
cal expressions for the Andreev bound states and supercuiTent, 
as a function of phase difference (p between the two supercon- 
ductors, and derived the 0-;7r transition by tuning the magnetic 
field. The transition takes place when the Zeeman energy Ez 
is of the order of the Thouless energy Ejh in the ballistic sys- 
tems. For the case of two conduction channels, we have ob- 
served a finite supercuiTent at ^ = (anomalous Josephson 
current) and direction-dependent critical current due to the in- 
terplay between the SO interaction and Zeeman effect. The 
critical cuiTent shows a cusp around the transition between 0- 
and TT-like-states, which is located at different positions for 
the positive and negative directions, as a function of magnetic 
field. 

Our model indicates the anomalous supercurrent in short 
junctions with more than one conduction channel, but not with 
single conduction channel. This is in contrast to the case of 
long junctions with L » ^, where the anomalous current is 
possible even with single channel."*^ However, we cannot ex- 
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elude that our result is specific to our model where the SO 
interaction works at a single scatterer 

Recently, the 0-;r transition and direction-dependent cusps 
of the critical current were observed in the Josephoson junc- 
tions of InSb nanowires when a parallel magnetic field is ap- 
plied.— A few conduction channels exist in the NWs, which 
is similar to the situation of our model, although L > ^ 
(L - 500 ~ 1000 nm, ^ ~ 350 nm) in the experiment and 
L <K ^ in our model. We are examining an extended model 
for L > ^ in which the scattering matrices, Sg and 5h, have a 
weak energy-dependence. Our preliminary result is not quali- 
tatively different from that presented in this paper concerning 
the case of two channels (anomalous current is possible with 
single channel in the model for L> In the experiment, the 
spin relaxation length by the SO interaction (^so ~ 200 nm) 
is comparable to the length of normal region L. However, it is 
hard to estimate the parameter pso in our model of two con- 
duction channels. We only know that pso - for ^so/L <K 1 
and 1 for f so /L» I. 
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In the case of two channels in the NW, the scattering ma- 
trix S in eq. ( fTSl ) is given randomly to follow the orthogonal 
ensemble in the absence of SO interaction and the symplec- 
tic ensemble in the strong limit of SO interaction. For the 
intermediate strength of SO interaction, the ensembles are in- 
terpolated with a parameter pso (0 < pso ^ 1), as described 
below. 



For the symplectic ensemble, the scattering matrix is writ- 
ten as a product of a diagonal matrix A and unitary matrix 
U, 



S = UAW 



(31) 



A is given by 



A = 



■^2" (8> 1 ; 



(32) 



with unit matrix 1 in the spinor space. Aj ( j - 1,2, ■■■ , 2N) 
are given randomly. The unitary matrix U is represented by 



(33) 



where 2N vectors ji^jj are complex. They are randomly cho- 
sen in such a way that iffj and f are orthgonal to each other 
for 1 <j,k< 2N. 

When a matrix in the symplectic ensemble is given, we 
make a matrix in the orthogonal ensemble as follows. From 
i^j, a real vector Xj is defined as Xj - Rei^j for spin compo- 
nent cr = + 1 and Xj - for spin component cr = - 1 . From 
jjCjj, an orthonormal set of 2N vectors |jEj| is created using the 
Gram-Schmidt orthonormalization. Then we obtain a matrix 
by eq. dSTI ) with U - (^Xi,gxi,- ■ ■ ,X2n,8X2n), where Xj and 
gXj have spin components of cr - +1 and -1 only, respec- 
tively. 

For the intermediate strength of SO interaction, we make 



^Xj+pso(l/fj-Xj). 



(34) 



We orthonormalize the vectors i^ifr'., gi^'*^ and construct U. 
This yields the scattering matrix in eq. (1311 1. 
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